Barkhausen-type noise in the resistance of antiferromagnetic Cr 
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PACS 73 . 50 . -h - Electronic transport phenomena in thin films 

PACS 75 . 30 . -m - Intrinsic properties of magnetically ordered materials 

PACS 75 . 30 . Ds - Spin waves 

Abstract -We present an experimental study of the changes generated on the electrical resistance 
R{T) of epitaxial Cr thin films by the transformation of quantized spin density wave domains 
as the temperature is changed. A characteristic resistance noise appears only within the same 
temperature region where a cooling- warming cycle in R{T) displays hysteretic behavior. We 
propose an analysis based on an analogy with the Barkhausen noise seen in ferromagnets. There 
fluctuations in the magnetization M{H) occur when the magnetic field H is swept. By mapping 
M ^ ^0 and H ^ where ^0 corresponds to the order parameter of the spin density wave, we 
generalize the Preisach model in terms of a random distribution of resistive hysterons to explain 
our results. These hysterons are related to distributions of quantized spin density wave domains 
with different sizes, local energies and number of nodes. 



Among simple metals, Cr is the only one that shows 
an antiferromagnetic electronic ground state correspond- 
ing to an incommensurate spin density wave (SDW). [1^ 
A single crystal of Cr cooled below the Neel temperature 
develops domains whose SDW wave vector Q may be ori- 
ented along any of the main crystallographic directions 
with the same probability. These domains are separated 
by walls which can move back and forth with temperature, 
or that can be suppressed by cooling in an applied mag- 
netic field H. [2] Using coherent x-ray diffraction Shpyrko 
et al. [3] measured the noise spectrum at different tem- 
peratures in Cr single crystals produced by domain wall 
motion separating domains with different Q orientations. 
Independently, Michel et al. [4] observed the presence of 
spontaneous fluctuations on the resistance of Cr films, re- 
sembling a telegraphic code. The resistance jumps were 
associated to the thermal motion of two different domain 
walls: rotations of Q domains, and rotation of the po- 
larization vector T] within a single Q domain. However, 
it has been shown that in thin films the SDW confine- 
ment leads to the quantization of the wave vector Qn, N 
being an integer, which orients perpendicular to the film 
surfaces. [5l[6] This quantization is due to the boundary 
conditions imposed by the film surfaces or interfaces. Rich 



hysteretic phenomena in the electronic transport proper- 
ties [71 [8] have been ascribed to the development of do- 
mains with different Qat corresponding to a different num- 
ber of nodes of the SDW. Consequently, in this case we 
would expect a different domains structure and dynamics 
compared to that seen in single crystals or thick films. 

In ferromagnetic materials, the existence of domain 
structures that can move due to the action of an exter- 
nal field H produce the well known Barkhausen noise |9] 
in the magnetization M{H). However, in antiferromag- 
nets or SDW systems the noise produced by domain wall 
motions as an external parameter is varied, has been more 
elusive to measure due to the lack of a net magnetic mo- 
ment. 

In this Letter we address an up to now unattended issue 
related to the resistance changes generated by the move- 
ment or transformation of antiferromagnetic domains that 
have different number of nodes (N or N+1) in the quan- 
tized SDW state present in thin Cr films. We show that 
the whole hysteretic behavior in R{T) is highly reminis- 
cent of that observed in ferromagnetic materials where 
the control field is and the measured magnitude is the 
magnetization M{H). In these systems, Barkhausen noise 
has been successfully explained in terms of magnetic hys- 
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terons p^QlfTT] using the Preisach model [TT] . In our trans- 
port measurements the control field is the temperature T 
and the measured magnitude is the film resistance R{T). 
Consequently we propose a mapping of H ^ T and of 
M ^ where is the order parameter associated 
with the SDW. Using this we show that a simple phe- 
nomenological extension of the Preisach model based on 
resistive hysterons [10] reproduce quite satisfactorily the 
experimental data. With this model we extract relevant 
conclusions about the origin of the hysteretic phenomena, 
the distribution of energy barriers and the processes in- 
volved in antiferromagnetic domain transformation with 
temperature. These findings may be relevant to appli- 
cations in spintronics, where the understanding of the dy- 
namics of antiferromagnetic domain walls is of importance 
to develop better pinning layers, as well as in general prob- 
lems of antiferromagnets in confined geometries. 

Cr thin films were grown epitaxially on MgO (100) sub- 
strates using DC magnetron sputtering. The films were 
characterized by x-ray diffraction showing rocking curves 
with an angular dispersion at FWHM around the [002] 
peak of 0.5° indicating a very good degree of epitaxy. Re- 
sults from AFM scans show a mean surface roughness of 
about 20A, which is the third part of the SDW wavelength 
at low temperatures. The samples were patterned in a 
four terminal configuration bar of 2 mm long and 140 /im 
wide using photolithography and chemical etching. The 
resistance was measured using a dc current of 10/iA as a 
function of T between 20 and 320 K in steps of 0.5 K in a 
commercial cryo-cooler. In each step the temperature was 
stabilized within 30 mK giving an intrinsic error in the 
resistance measurement between 0.3 to 0.8 m^l, almost in- 
dependent of T. Each R{T) curve took about 24 hours of 
measurement. 

In Fig(T]we present our main experimental findings. The 
changes in the resistance produced by the cooling and 
warming of the sample are small jumps (|Ai?| ~ 10~^ ft). 
Since the resistance is of the order of tens of Ohms (i?(300 
K) 63 r^) a plot of R{T) shows no appreciable features. 
Instead, in panel (a) we show an amplification of the effect 
by plotting the difference of the resistance derivatives for 
a complete (20 K to 320 K) cooling (c) and warming (w) 
cycle: dRc/dT - dR^/dT for a 550 A thick film. The solid 
line corresponds to a smoothing of the raw data using the 
Savistky-Golay method. The inset shows the difference 
between the raw and smoothed data. Although noise is 
seen in the whole temperature range, it is evident that its 
amplitude strongly increases in the temperature window 
where the irreversible behavior in R{T) sets in (see inset). 
This clearly indicates that the observed noise in dR/dT is 
not exclusively due to thermal fiuct nations. We claim that 
it has to be related to the same mechanism that governs 
the hysteretic behavior. 

Neutron diffraction [12], x-ray diffraction [7,, and elec- 
trical transport measurements [7l[8], show that the hys- 
teretic behavior in R{T) and dR{T)/dT is due to the exis- 
tence of antiferromagnetic domains with different number 



of nodes N in the confined SDW. These domains switch 
from to + 1 as the temperature is decreased. As 
we will show and elaborate in more detail below, we as- 
sume that the noise pattern seen in Fig(T] (a) is mainly 
due to the switching and growth of these domains as the 
temperature of the film is changed. From previous re- 
sults [8] we know that the hysteretic region enlarges in T 
as the film thickness increases. Taking this into account, 
we performed measurements of cooling- warming cycles of 
different amplitudes in T in a Cr film with a thickness 
of 11 00 A for which the irreversible behavior embraces a 
wider temperature range as compared to that found in 
thinner films. The cycles were measured starting always 
at a temperature of 320 K, cooling down the film to a given 
final temperature (which we call Tj), and then warming 
it up again to 320 K. The results of such measurements 
plotted as AR{T) = Rc — Rw are shown in panel (b) of 
FigHJ Clearly, as Tf is increased the cycles reduce their 
amplitude, the maximum shifts to higher temperatures, 
and the curves tend to coincide with the complete cycle 
curve at high enough temperatures. From these data we 
may construct the equivalent to the first order reversal 
curves (FORCs) [13] in ferromagnets. They are shown in 
panel (c) and are defined as ARs = Rc — R1/2 cool- 
ing cycle and ARs = Rw — R\/2 warming cycle, 
where ^^^72 ~ (^c + ^2;)/^ is taken from the complete 
cycle {Tf = 50 K). 

The measured resistance depends on the domain distri- 
bution, the domain wall structure and the energy barriers 
in a complex free energy landscape. Due to the problem 
complexity we resort to a simple phenomenological model 
to describe the physics behind the observed hysteretic be- 
havior. We assume that the building blocks of the model 
are resistive hysterons that correspond to small hystere- 
sis cycles in R{T). These are the equivalent of magnetic 
hysterons commonly used to explain M{H) hysteresis in 
ferromagnets. There, the Preisach model naturally in- 
corporates the fact that the magnetization increases with 
magnetic field and uses the magnetic hysterons to account 
for the increase (decrease) of the total magnetic moment 
as the field increases (decreases). In our case, as we show 
below, the resistive hysterons should describe the tran- 
sition from a high to a low resistance state as the tem- 
perature is decreased in order to properly account for the 
experimental data. From the detailed analysis of Kumma- 
muru's experimental data on electrical transport [7] it was 
inferred that the state with N nodes has larger resistance 
(or a smaller number of effective carriers) than the state 
with N+1 nodes. Our recent measurements of the Hall 
coefficient as a function of temperature also support this 
view showing hysteretic behavior with a higher number of 
carriers while heating [M]. Although this seems counter- 
intuitive (in the low temperature phase we would expect 
the amplitude of the SDW to increase and consequently to 
reduce the number of carriers), numerical simulations (in 
multi-layer systems with magnetic interfaces) show that 
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Fig. 1: (color online) Experimental data: (a) Noise in the difference of derivatives dRc/dT—dR^/dT associated with the complete 
loop in a 550A thick film. The solid line corresponds to a smoothing of the data (see text). The inset shows the difference 
of the data and the smoothing, (b) Difference AR(T) of the resistance in cooling-warming cycles of different amplitudes in T 
(with final temperature T/), for a llOOA thick film, (c) FORCs, LO (T/ = 50 K), LI (T/ = 150 K), L2 (T/ = 185 K) and L3 
(Tf = 225 K), calculated from the data in (b) with ARs as defined in the text. 



the amplitude of the order parameter is drastically reduced 
when the number of nodes increases by one [15j. Based on 
this evidence, we consider that domains with N+1 nodes 
have lower resistance than those with N nodes. With our 
phenomenological description we attempt to account for 
the fact that the resistance is given by the presence of 
domain walls and mostly by a contribution from the do- 
mains. Hysterons represent a transition from states with 
N (high T) to N+1 (low T) nodes. A second source of re- 
sistance noise could be the redistribution of domain walls 
as growth of the domains implies loss of domain- wall scat- 
tering. The shape of the resistance hysteresis loop for the 
complete cooling- warming cycle shows that the dominant 
effect is the change in the domains resistance: Rc is always 
larger than R^ . If the irreversibility were to be dominated 
by the domain- wall scattering of pinned walls that jump as 
the domains grow, the shape of the resistance loop would 
be different and we should expect a temperature interval 
where R^ > Rn. 




^ T[Z.u.] 



Fig. 2: (color online) Theoretical model: (a) difference of the 
resistance in cooling-warming cycles of different amplitudes in 
T. Inset: resistive hysteron. (b) FORC curves LOM, LIM, 
L2M, and L3M obtained from the data in (a) with ARs as 
defined in the text. 



R(T) = aT^b- 



tanh 



T-To^± D, 



a 



(1) 



We characterize these hysterons with three independent 
parameters: the center of the loop Tq, the half- width of 
the loop D and the amplitude of the change in resistance 
Z. See inset in Figl2](a). 

Since we are interested in the derivative of R(T) 
it is better to mathematically describe the hysterons 
using an analytical function. We choose r(T) = 
Z [tanh( ^~^Q^^ ) + l] , with a < D, that describes the 
switching between 2Z and at Tq =F In order to take 
into account the complex domain structure, we use a col- 
lection of independent hysterons with random parameters 
{Tq^D^Z) each. Variations stand for different domain 
sizes, local energies involved and different number of nodes 
(N or N+1). Once defined a collection of Nh hysterons, 
the total resistance at a given temperature R{T) is just 
the sum of the resistance of each hysteron ri{T). We can 
then write 



where the + (— ) sign stands for cooling (warming), and 
we have added a linear trend aT + 6 as the base where 
hysterons are mounted. This is a reasonable assumption 
to account for the phonon and impurity scattering con- 
tributions to the resistance in the temperature range of 
interest. 

The results obtained from this model are shown in Figl2l 
Based on the shape of the FORCs we adopted a random 
gaussian distribution for each of the parameters Tq^D^Z 
with average values and standard deviations: (Tq) = 7, 
cTTo = 0.6; {D) = 0.5, ao = 0.12; and (Z) = 10"^, 
az = 0.15 X 10~^. In panel (a) cooling- warming cycles 
have been computed using different distributions of hys- 
terons for each cycle. There are two main effects of on 
R{T), namely the enhancement of the difference between 
Rc and i?^ in the hysteresis region, and the smoothing 
of the curve shape. When is very small, ~ 10^, 
discretization plays a mayor role, giving a steps-shaped 
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curve. To reproduce the features of our experimental data 
we took Nh 700. In Panel (b) we show the FORCs 
obtained with the model. 

A comparison of Figs|T] and [2] nicely shows that the re- 
sults obtained with the model are in very good agreement 
with the experimental data. This supports our view of the 
hysteresis process as a complicated movement of domain- 
walls with no memory of the path followed to reach the 
state with N + 1 nodes by lowering the temperature, or 
with N nodes by increasing it. 

In Figl3]we display the noise in the loops LO to L3, that 
is the peak structure present in dAR/dT. Panel (a) shows 
the experimental data shifted vertically by an arbitrary 
offset for the sake of clarity. In panel (b) the noise pat- 
terns computed using the model are plotted. A small ran- 
dom gaussian noise has been included to account for the 
experimental thermal fluctuations in R{T). The similarity 
is remarkable and we may say, using the analogy with fer- 
romagnets, that this noise in R{T) is the resistive version 
in SDW systems of the Barkhausen noise characteristic 
of ferromagnetic materials. [TBI However, we expect some 
differences between the ferromagnetic and SDW systems. 
In particular, the long-range dipolar interaction present in 
ferromagnetic films plays an essential role in the domain 
dynamics. Unfortunately, we do not have a direct access 
to observe the dynamics. Since in our case the control pa- 
rameter is the temperature, the dynamics of the domain 
redistribution is affected by the way the sample thermal- 
izes after each temperature change. Therefore, we only 
have access to the final state at each temperature. 

Figsl3] (a) and (b), seem to indicate that each experi- 
mental noise pattern is different. 

In order to give a more qualitative analysis of this state- 
ment we digitalize the data and introduce a function that 
allows us to obtain a code-bar type identification of each 
pattern. We create our code-bar in the temperature axis 
by dividing it in M small intervals T^. We subtract from 
each noise pattern in Fig. 3(a) the smoothed curve to ob- 
tain a pattern as the one shown in the inset of Fig. 1. 
Taking the highest (positive or negative) peak amplitude 
in the resulting pattern we define the relevant scale that 
is divided in intervals Wk of fixed size A. The interval Wk 
extends from kA to {k + 1)A with A; = 0, 1, N. For the 
pattern of each loop Li and each window Wk, the code 
function is defined as: 

J Codei^k[Tn] = 1, if the peak amplitude G Wk 
\ Codei k[Tn] = 0, otherwise 

(2) 

In this way, for each pattern and window we have an 
M- dimensional vector whose components are or 1. The 
square modulus G^^ = ICodei^klTnW^ = Codei^k[Tn] • 
Codei^k[Tn] of each code gives the total number of peaks 
with amplitude in the window Wk. In order to compare 
the noise spectrum of the different loops we first note that 
although all codes have the same dimension M, the first Ni 



components of loop Li (i=l,2,3) are zero by definition. In 
Fig. [31(c) we show the normalized functions G^^/ (M — Ni) 
(upper curves) that are equivalent to the histograms of 
the noise amplitude. The solid line is a fit with a Gaus- 
sian function. In order to analyze the reproducibility of 
the noise in different loops (Figl3](a)), we define the corre- 
lation function as the scalar product of two code functions, 
i.e. G^j = Codei^k[Tn]-Codej^k[Tn]/{M-Nj) with i < j. 
The lower curves correspond to Gqi (full circles), G02 (full 
diamonds) and Gqs (full squares). The smaller correlation 
seen in all Goi correspond to stochastic correlations among 
two different gaussian distributions, one for loop LO, and 
the others for loops L^, with i = 1,2,3. Correspondingly, 
the solid line is a fit with the same upper Gaussian func- 
tion but squared. 

It is important to realize that if the noise patterns were 
exactly the same for all loops, irrespective of its tempera- 
ture span, all curves should coalesce onto a single curve (as 
clearly is the case for the autocorrelation functions Ga). If 
instead, only part of the loops were correlated (for exam- 
ple the common path when cooling down) the G^j should 
give larger values than the correlations observed among 
two different gaussian distributions. Therefore, this result 
illustrates in a more quantitative way that each realization 
of the experiment follows a different path with respect to 
the distribution of SDW domains, giving rise to differ- 
ent resistive noise patterns in each different run. Finally, 
we mention that the noise amplitude with quasi- Gaussian 
distribution is present for all studied films, and has a stan- 
dard deviation that increases with film thickness. 

Another way to analyze in more detail the noise pat- 
terns, would be to construct a FORC diagram [17j with 
the experimental data. This procedure would allow to 
estimate the distribution of hysteron parameters directly 
from the experiments. However, this task would demand 
more than a year for measuring a reasonable number of 
loops, say 200, because each loop takes almost 50 hours of 
measurement. The good agreement between experiments 
and simulations led us to evaluate the FORC using our 
simple model. The resulting FORC diagram is similar to 
the one obtained for systems of single domain ferromag- 
netic particles (SDFP) with a narrow size distribution. In 
these systems each particle is subject to the effect of and 
external magnetic field and to a random field Hint due 
to the dipolar interactions. In view of these similarities, 
we can make an analogy between the Cr SDW resistiv- 
ity hysteretic behavior and the much simpler case of the 
SDFP. In the SDFP system each hysteron corresponds to 
the orientation of one particle magnetization along the ex- 
ternal field direction. The random fields play an important 
role in determining the shape of the FORC diagram. In 
the case under study here, hysterons correspond to tran- 
sitions between different configurations of domains with 
N or N+1 SDW nodes. The role of the random fields is 
played by a random distribution of Tq which is crucial in 
order to reproduce the experimental observations. Finally, 
the particle size distribution in the SDFP case is here rep- 
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Fig. 3: (color online) Noise in dAR/dT for the loops shown in: (a) Fig[T] (experiment), and (b) Fig|2] (model). Panel (c) shows 
the correlation function Gij = G{Codei, Codej) /min{Ni, Nj} as a function of the window position. See text 



resented by the distribution of resistance jumps Z, and 
the values of D are related to the energy barriers strength 
between different domain configurations. 

In summary, we have shown experimentally that the 
resistance of thin Cr films displays a characteristic noise 
in the same temperature region where hysteretic behav- 
ior is seen. We presented an extension of the Preisach 
model and the introduction of resistive hysterons in a sim- 
ple phenomenological model which reproduced quite sat- 
isfactorily the noise patterns in dR/dT and hysteresis in 
R{T). These are mainly generated by the switching of 
a random distribution of SDW domains with N or N-|-l 
nodes in the spin density wave. The main assumption 
supported by previous experiments [7] is that the domains 
with N+1 nodes have a lower resistance than those with N 
nodes. Our results present new evidence supporting this 
scenario. However, to our knowledge a microscopic theory 
is still lacking and new theoretical efforts are needed to 
complete the understanding of this interesting issue. We 
show that each realization of cooling-warming cycles fol- 
lows a different path with a different random distribution 
of domains. We conclude that the system presents a nar- 
row distribution of Tq that suggests a narrow distribution 
of met a stable domain configurations. Our approach also 
opens the possibility to study domain evolution as a func- 
tion of temperature in other antiferromagnetic materials 
in confined geometries. 
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